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SUHHABY 


The  two  dloMnelonel  problem  of  the  interaction  of  e  plene  week 
shock  wsve  with  e  cyllndrlcel  obstecle  of  erbltrery  cross  section  is 
considered.  An  Integrel  equation  for  the  surface  values  of  the  potwi- 
tlel  is  fomuleted  and  solved  epproxlmetely  for  the  case  of  e  square 
box  with  completely  rigid  boundaries. 


XHTRODUCTION 


This  report  la  coaearoad  vitli  cha  too  dlMooloMl  probloa  of  tha 
detaralnation  of  tha  praaaora  and  oalocity  fialda  raaoltiog  froa  tbo  lii« 
taractlon  of  a  plaoa  oaak  ahock  waoa  with  a  eyllodrleal  ohataela  of  arhl« 
trary  eroaa  aactlon.  toeh  problaao  haoa  baan  pravlooaly  conaldarad  and 
analytical  solutlona  for  a  naidMir  of  alaple  apaclal  gaoaatrlaa  hava  baan 
found  by  varloua  taebnlgnoo. 

Tha  aathod  of  laagaa  baa  baan  eoployad  to  conatrnet  a  aolutlon  for 
an  Infinite  plana  obataela  (1).  Separation  of  oarlablaa  baa  been  need  for 
•neb  alapla  gaoaMtrlea  aa  tha  circular  cylinder  (2).  lha  problao  of  an  ln> 
finite  wedge  baa  baan  aolwad  by  conforoal  Mapping  with  the  uaa  of  jaoMatflcal 
acouatlca  to  dateraina  tha  proper  region  of  aolutlon  (3).  This  latter  pro- 
blaa  has  also  been  solved  for  the  sMra  general  case  of  curved  wave  fronts 
by  tha  uaa  of  tha  appropriate  Green's  functions  (k).  Configurations  Invol¬ 
ving  sharp  corners  are  character  load  by  tha  fact  that  these  corners  act  as 
canters  of  diffraction  leading  to  different  solutions  In  different  regions 
of  space-tine.  This  feature  allows  the  construction  of  the  solution  for  a 
box  froa  that  for  an  Infinite  wedge  but  tha  procedure  Is  lapractlcal  for 
SMra  than  two  corners. 

Apart  fron  such  siapla  gaonetrias  tha  general  problaa  does  net  appear 
aaenabla  to  analytical  solution.  The  use  of  nuaarlcal  nathods  based  on 
finite  differences  to  solve  the  differential  equations  Is  lapractlcal,  evan 
with  large  scale  coaputers,  because  of  the  three  dlannsional  (two  space,  one 
tlaa)  character  of  tha  problaa  and  tha  presence  of  discontinuities  In  tha 
field. 

A  significant  slaplificatlon  can  be  aade,  however.  If  tha  problaa  Is 
rafonsulatad  In  tarns  of  values  of  tha  psaasMiia  on  the  surface  of  tha  ob¬ 
stacle  alone.  The  spatial  diaanslons  of  tha  problaa  are  reduced  froa  two 
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to  one  and  thla  roaalnlng  apace  dlMMalon,  If  the  obataela  is  finite,  ia 

t 

llaited  in  extent.  This  foraulation  is  particularly  suitable  for  physical 
problesM  uhich  require  only  a  solution  on  tbs  surface  of  the  obatacle.  If 
the  aolution  in  the  field  is  also  desired  it  can  be  constructed  ia  a  atraipht 
forward  aannar  froa  the  surface  values. 

The  procedure  developed  here  forai>lates  en  integral  equation  for  the 
the  fsisSdre’  et  an  arbitrary  field  point  in  teras  of  the  initial  wave  and 
an  integral  of  tiae  retarded  values  of  the  pneasefe  end  its  derivatives 
over  the  surface  of  the  obatacle.  An  integral  equation  on  surface  valuea 
alone  is  obtained  by  allowing  the  field  point  to  approach  the  surface  of 
the  obstacle;  valuea  throughout  the  field  aay  be  obtained  by  direct  inte¬ 
gration  over  the  surface  valuea. 

The  effect  of  the  discontinuity  in  the  pressure  and  the  velocity  at 
the  wave  front  on  the  aurface  can  be  separated  froa  the  resaining  surface 
effects  and  integrated  directly.  The  raaainlng  integrations  are  approxiastad 
by  assuaing  the  surface  psassura  ,to  have  an  average  value  over  specified 
steps  in  space  and  tiae.  The  intagrationa  are  Chan  replaced  by  suanations 
which,  because  of  the  tiaa-retarded  effect  lead  to  succassiva  algebraic  non- 
siaultanaous  aquations  on  the  unknown  surface  values.  The  fact  that  these 
aquations  are  not  siaulggnaaos  i*  vital  because  it  peraits  the  use  of  a 
large  nuabar  of  assh  points  without  prohibitive  coaputations. 

The  pressure  distribution  on  tbe  surface  of  a  rigid  square  box  undar 
a  syaaatric  plana  pulse  loading  is  found  by  desk  coaputations  for  tiae 
steps  up  to  one  transit  tiaa.  The  portions  of  the  solution  corresponding 
to  the  infinite  wedge  give  excellent  agraaaant  with  the  known  analytical 
solutions  (3). 
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1,  Initial  -  Boundary  »  Probli  FemulatlOB. 

The  praaaura  flald  that  raaulta  froa  tha  Intaraetion  of  a  plana 
acouatlc  shock  with  a  two-diasnsional  obstacla  of  arbltnayduipOy  rigid* 
and  fixed  in  an  acoustic  ■adiua  is  to  ba  dstahsinad  f or  tha  purpoaa  of 
snalpals  it  in  aora  convaniant  to  consider  a  finite  cylindrical  (Ostade  of 
arbitrary  length  and  of  constant  cross  section  rather  than  a  two  diaansional 
obstacle.  The  field  corresponding  to  this  three  diaansional  problca  is  iden¬ 
tical  in  every  plane  section  noraal  to  the  axis  of  tha  cylinder,  i.a.  inde¬ 
pendent  of  distance  along  the  axis  of  the  cylinder,  for  tines  less  than  th^ 

I 

propagation  tine  froa  tha  ends  of  the  cylinder  to  the  nearest  section  con¬ 
sidered.  Over  these  tiae  Intervals  the  three  dlnensional  field  provides 
the  desired  two  diaansional  solution.  This  is  a  consequence  of  the  hyper¬ 
bolic  character  of  the  wave  equation  which  governs  tha  field. 

In  this  three  diaenslonal  field,  let  a  cartesian  coordinate  systea 
x',  y*,  s'  be  introduced  with  the  s'  axis  taken  parallel  to  tha  generators 
of  the  cylinder.  The  cross  section  s'  ■>  constant  is  bounded  by  an  arbitrary 
piece-wise  saooth  curve  x'(s),  y'(s)  where  the  paraaater  s  is  arc  length 
along  the  boundary  of  the  cross  section  [fig.  1  ].  At  t  ■  0  tha  obstacla 
is  subjected  to  a  plane  acoustic  shock;  the  line  of  contact  between  the 
shock  and  the  obstacle  is  taken  as  the  s'  axis.  Tha  Incident  shock 
raprasants  a  step  in  pressure  whose  aagnitude  is  used  to  noraulise  tha 
pressures  idilch  are  given  as  differences  froai  tha  prassura  in  tha  imdisturbad 
atata  ahead  of  tha  shock.  Consequently,  the  stataa  ahead  and  behind  tha 
shock  at  t  a  0,  correspond  respsetively  to  p  *  0  and  pal  thus  defining 
tha  initial  stats. 

♦rtie  Mtiiod  developed  In  this  report  is  applicable  to  the  smre  general  case 
of  a  non-rigid  obatacla;  see  fortheoaing  reports. 
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Thu  interaction  field  nay  be  described  In  terns  of  a  continuous 
velocity  potential  which  satisfies  the  wave  equation  and  is 


related  to  the  pressure  through  p  a  -f  o 


dent  shock  is  defined  by 


f.  ifCwAVt)  5 


In  terms  of  the  inci> 


o 


;  xcos^t';4^P* 


where 


X*  ^/c  ,  ^  =  ,  i-  z'/c. 


^ being  the  speed  of  sound  and  ^  is  the  angle  between  the  x  axis  and  the 

normal  to  the  incident  shock  front.  This  together  with  the  rigid  surface 
(D 

condition  — — •  0  ,  n  is  the  normal  to  the  obstacle  surface,  estab¬ 
lishes  the  following  initial-boundary  value  problem  for  . 


D.t. 

I.C. 


b.c 


a'lp  =  «PkK  +  >i>yy  ^  =0 


n  n  (  t-xa>sd -jbin/S  ;  o> 

Jo  yKy,i,o)=  -j 

O 


s 


m 

I  \ 

I  0 


(l.l) 
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Although  equations  (1.1)  spoclfy  the  probloa  coaplotaly,  the  fora 
is  not  aultablo  for  ahalysio  nhan  dMllng  with  arbitrary  ahapas.  InatMd 
an  oquivalant  foraulation  in  taras  of  an  integral  aquation  can  be  developed 
which  is  sore  aaenable  to  further  treatnent. 


2.  Integral  Equation  Porraletion. 

The  integral  equation  can  be  developed  by  treating  the  problea  as  a 
characterlatic«boundary«value  problea  in  It^'diaenslonal  space-tiae  (x,y,s,t } 
to  which  Green's  identity  aay  be  applied,  rather  than  as  an  initial>bound« 
ary-value  problea.  Consider  the  surface  foraed  by  the  intersection  of  the 
characteristic  hyper-plane  representing  the  incident  frOatiand  the  cylin¬ 
drical  hyper-surface  representing  the  obstacle,  in  space-tiae.  The  union 
of  the  Influence  doaains  of  all  the  points  of  the  surface  of  intersection 
is  a  doaain  in  h-space  exterior  to  which  the  "secondary"  (disturbance)  po¬ 


tential  defined  by  <^s  ■  ^w  vanish.  The  boundary  of  this 

region  consists  of  reflection  and  diffraction  characteristics.  For 
to  be  a  continuous  solution  for  all  (x,y,s,t),  aust  also  vanish  on 

this  boundary. 

In  addition  to  the  characteristic  surface,  say  ^  (’‘o*^o'*o*^o^' 
corresponding  to  the  secondary  potential,  there  exists  at  any  point  (x, 
y,s,t),  a  characteristic  half-cone  T*  : 

t-to  *  m  R 


These  characteristic  surface,  together  with  the  obstacle  surface  fora  a 
closed  regionVin  h-space  idiose  projection  in  physical  space  is  exterior 


to  the  obstacle 


-6. 


For  the  nave  operator  Q  S  +■  ^  JJl  ““ 


Green 'a  Identity  tekea,  the  form 


where  repreaente  conorael  differentiation.  If  now  rT  la  identified 

t-t 

with  ^  g>  u  taken  aa  the  progreaaive  wave  aolution  u  ■  1-  -■”  ,  and  V,  S 
correspond  to  the  above  nentioned  region  and  bounding  surface,  the  identity 
reduces  to 

c* 

C  £  corresponds  to  a  hyper  •cylinder  of  radius  ^  cut  out  of  the  region 
V  to  account  for  the  singularity  in  u  at  r  ■  0;  S'  denotes  the  surface 
of  the  obstacle  plus  the  surface  of  the  secondary  disturbance  in  k-space. 
For  the  chosen  function  u,  the  integrals  over  C  ^  can  be  reduced  to 


( trX,o)  dlbo 


In  addition,  on  that  portion  of  S'  corresponding  to  the  secondary  front, 

^  ^  >  0  and  consequently  *  0  since  the  front  ,  is  characteristic. 

On  the  obstacle  surface  .  Hence,  eq.  (2.1)  sian 


plifies  to 


=  5^  Hi  (a 

i  *  $* 


Differentiating  both  sides  twice  with  respect  to  t,  pidlds 

/  i- 
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Here  the  U«space  surface  Integration  has  been  separated  Into  a  surface  In¬ 
tegration  S  over  the  physical  surface  of  the  obstacle  with  linlts  de- 
0 

tensined  by  the  intersection  of  the  cone  J"*  and  the  secondary  front  on 
the  surface  of  the  obstacle  in  h-space  and  a  time  Integration  whose  upper 

limit  Is  the  cone  r. 

The  tine  differentiation  may  be  carried  through  the  S  integration 

o 

since  the  integrand  vanishes  along  the  curves  of  the  intersection;  u  vanishes 
along  i  while  ^  vanishes  along  the  secondary  front.  Hence  ^ 

Since  the  solution  function  ^  ^  is  uniquely  determined  by  the  specifica¬ 
tion  of  alone,  as  given  by  eq.  (1.1),  it  follows  that  eq. 

(2.2)  represents  an  integral  equation  for 

A  significant  simplifications  can  be  achieved  by  establishing  that 

))  I  «  R*  ^'^5^0  R  ^  (2.5) 

S0  ia*t>R 

2.^  3 

for  then  eq.  (2.2)  can  be  written  in  terms  of  the  total  potential  z  * 


1.  This  can  be  considered  a  generalization  of  the  classical  Klrchoff's 
formula,  the  latter  dealing  with  a  physical  space  of  fixed  extent, 
l.e.  independent  of  time. 

2.  An  attempt  to  derive  this  eq.  for  directly,  would  involve  consider¬ 
ation  of  integrals  over  the  characteristic  secondary  front. 

3.  An  alternate  heuristic  derivation  of  this  equation  by  means  of  Dirac 
delta  functions  is  given  in  Appendix  I. 
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Es.  (2.3)  say  be  verified  by  extending  the  definition  of  ^  or 

w  into  a  space  time  region  whose  projection  in  physiesl  specs  is  interior 
^  s 

to  the  obstacle,  by  means  of  a  saltus  problem.  Consider  the  region  9  in 
4-spaee  bounded  by  the  three  characteristic  surfaces  correspohding  to  the 
secondary  front,  the  cone  T*  and  the  incident  wave  front.  This  region  con* 
tains  an  interior  boundary,  namely  the  surface  of  the  obstacle.  For  V  con¬ 
sidered  as  a  single  region  it  is  possible  to  define  a  solution  ^  which  is 
discontinuous  but  tdiose  normal  derivative  is  continuous  across  the  obstacle 
surface. 


This  solution  is  identical  to  ^  in  V  as  defined  by  eq.  (2.2)  while 

it  vanishes  in  the  region  V-V.  That  the  latter  is  the  correct  extension  follows 

from  the  condition  ^  s  0  on  the  incident  front  and  ■  0  on  the 

obstacle-surface.  This  in  turn  implies  fV  ■  •  ^  in  %V.  Consequmtly 

(i  ^*^0 

Green's  Identity  loay  be  applied  to  this  region  with/ir«  j  end  n  •  1  -  ■■■  "y". 
Since  Che  field  point  x,y,z,c  is  exterior  to  V-V,  eq.  (2.?)  follows. 


Returning  to  eq.  (2.4),  let  it  be  differentiated  with  respect  to 


t  and  multiplied  by  with 


>9  /  >  n  set  equal  to  sero: 


I* 


5) 


where 


and  '^oi*  ,  5*R  are  the  roots  of  (^sh*^)  “  *  O 

The  differentiation  may  be  carried  out  to  yield 
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f) 


♦ 


47r  ^ 

s,t  ^  4»*t*ie 


(2.6) 

) 

&•  &»>E 


Here  R^  ■  (x-x^(8)))^  +  (y-y^C*))^  +  Z^  ,  so  that  p  is  independent  of  s. 
Ublike  the  situation  encountered  in  obtaining  eq.  (2.3),  the  integrands  of 
eq.(  2.5)  do  not  vanish  at  m  -  z.  This  gives  rise  to  the  last  tern 

in  eq.  (  2.6)  %)hich  involves  values  of  and  p  ianediately  bAind  the 

wave  front  (1^  .  '**  interpreted  as  the  effect  on  a  field  point 

of  the  discontinuity  in  the  incident  pressure  field;  it  can  be  evaluated 
readily  for  a  given  x(s),  y(s)  since  ^  -  0  and  p(x^-s)  is  known. 

Eq.  (2.6)  can  be  interpreted  as  an  integral  equation  for  the  pres¬ 


sure.  Consider  its  behavior  as  the  field  point  (x,y)  approaches  the  surface 
of  the  obstacle.  The  integrands  become  infinite  at  the  point  R  a  0.  By 
considering  the  behavior  of  the  integrals  in  the  vicinity  of  this  point,  it 
is  readily  established  (appendix)  that  these  iaproper  integrals  converge  in 
the  ordinary  sense,  in  fact,  the  nnkxlbution  S»  the  double  integral  from  the 
inediate  vicinity  of  this  point  approaches  the  limit  value  1/2  p(x(s),  y(s),t). 
The  remaining  portion^  from  the  double  integral  is  finite  and  in  the  special 
*•••  curface  point  located  directly  behind  the  wave  front,  this  contri¬ 

bution  vanishes.  In  this  case  the  single  intstral  alt*  vanishes.  This  is 
in  accord  with  the  well-known  results  that  for  the  interaction  betwaaa  a  plana 
wave  and  an  arbitrary  obstacle,  the  pressure  imnsdiately  behind  the  refloetad 
wave  front  on  the  surface  of  the  obstecle,  is  Just  twice  the  pressure  bdhind 
the  incident  wave. 
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Eq.  (2.6)  may  than  ba  wrlttan 

lha  proeadura  aiiq>loyad  in  solving  aq.  (2.7)  for  p(s,t)  involvaa  tha 

approximation  of  tha  double  integration  as  follows.  Let  p(8,t)  be  assuoMd 

to  have  a  constant  mean  value  over  fixed  intervals  in  a  and  t.  Associated 

with  these  intervals  are  specific  regions  in  s^  and  z  in  which  p  is  constant. 
Bo  ** 

In  addition  let  ^  «  be  replaced  by  a  backward  difference  in  Cine.  Then  tha 
®  o 

double  Integral  can  be  approximated  by  a  double  summation  of  time  retarded 

values  of  p,  whose  coefficients  are  integrals  of  4  "V"  ^  and  -4 

over  tha  associated  regions  in  s  and  z  .  These  latter  integrals  can  ba 

o  o 

evaluated  for  a  given  obstacle  shape.  With  this  approximation,  eq.  (2.7) 
can  be  applied  to  successive  intervals  in  s,  t  to  yiald  a  set  of  successive 
algebraic  equations  for  the  mean  values  of  p.  The  details  of  this  SMthod 
are  illustrated  in  section  ^  where  a  box>shaped  obstacle  is  considered. 

Computations  have  indicated  that,  since  the  major  contribution  to  p  is 
made  by  the  terns  in  eq.  (2.7)  which  can  be  evaluated  exactly,  the  values  of 
P(»#t)  obtained  in  this  manner  are  relatively  insensitive  to  the  magnitude 
of  the  fixed  intervals  (provided  ^  !)•  It  is  anticipated  that  this 
will  be  the  case  in  general. 


3.  Solution  for  the  Case  of  a  Box-Shaped  Obstacle 

The  specific  configuration  for  which  eq.  (2.7)  has  been  solved  Is  a 
square  box  whose  boundaries  are  defined  by  S^:  x  ■  0,  S^t  y  >  0,  S^t  y  ■  a, 
S|, :  X  ■  a.  The  direction  of  notion  of  the  incident  wave  is  chosen  for 
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eonvenlmc*  of  iTMOtry,  to  bo  ouch  that  tha  noroal  to  tha  front  eoineidaa 
with  a  diagonal  of  tha  aquara,  l.(  3  r  .Jk  (fig.  Z  }.  In  thin  eaaa  aq. 


(2.7)  takaa  tha  fora 


f  u,t)  r  a  ^  *  i  ^  I  j/*-  [ 

R*  (fiL-a  ('**'*■  r  Sfes-  JL  1 

with  >ma5Jw»«a» 

The  aiBbol  i  rafara  to  tha  aurfaca  containing  tha  point  a,  at  idiieh 
tha  preaaura  ia  p^.  ia  not  included  in  tha  intagrationa  ainca  tha  ox- 

preasiona  containing  vaniah  for  a  aurfaca  of  conatant  curvature  (ap- 

pandixJI  ).  In  tha  axp^aaion  for  ,  tha  tana  p'^  ia  equal  to  tha  con¬ 
atant  two  for  J  >  1,  2  and  ia  equal  to  aero  for  J  >  3,  ^  (aaa  appendix  IL  ). 


The  regiona  of  integration  in  a^,  a^  on  each  aurfaca  are  dataminad  by  tha 


elllpaea: 

5.. 1  :  t-  “/JI  ■  "R 

5., H  :  t  -  W  -  Vc  ‘  "R  (H  X,d) 


with  the  raatructlona  that  for  a  >  a,  a  ia  eat  equal  to  a  and  for 
a^j^  ^  Of  a^^  ia  aet  equal  to  aero. 

The  double  integral  ia  approxinated  by  aaauaing  p'^(',t}  to  have  the 
conatant  value  pj^  (t-j(-fl)  over  the  apace-tiae  interval  (k-l)/2  ^a^^k  /S  , 
t  -j(  <t  ^  t  -  5^  -i-l).  The  time  derivative  ^ ia  approxinated  by  a  back- 
warda  difference  in  tine,  pj^  (t-l-t-l)-  pj^  (t  -ft)*  for  pointa  a^  awapt 
over  by  the  incident  front  within  the  tine  interval  being  conaldmred,  the  da- 
rivative  ia  replaced  by  p|[(t^)-2  for  J  ■  1,  2  or  for  J  ii  3,h,  Tha 

integration  nay  than  be  replaced  by  a  double  auanation  over  k  and  X  in 
which  tha  coafficianta  of  p^  (t^)  are  known  intagrala. 

It  ia  alao  convenient  to  nodify  alightly  tha  regiona  of  integration  in 
t^f  in  order  to  reduce  tha  nunbar  of  coafficianta  of  p^  that  need  to 
ba  conputad.  For  thia  purpoaa,  tha  regiona  of  int^ration  are  approxinated  byt 


i-', 
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vlth  th«  prtvlous  rMtrletiOQS  on  ttill  valid  and  balng  daflMd 

aa  aaro  in  tha  raglon  axtarlor  to  tha  allipaaa  t  -  a^//?  ■  t,  t  -  a/V?  * 
l.a.  ahaad  of  tha  mmm  front,  lha  eoaffieianta  obtalnad  In  thla  aaimor  ara 

In  fact  tha  aaaa  for  any  orlantation  of  tha  Ineldant  front  ralativa  to  tha 

box. 

It  than  follotia  that  tha  approxlaata  ayatam  to  ba  aolvad  la  glvan  byt 


4  t  P 


f>)*  ill 

«  <1!  k->  I  *  ‘  “ 

rx - r - .  f  (2.7) 


i  i;-  fe}. 


j)  d  •«  d 

|tt  paaaad 

ovar  point  k  batvaan  t-|  and  t-jKfl;  othanriaa 

fk*  (t-Avi)  —  X  oA 

f»  Ci-t+O  -  O 


tfhere  P  la  tha  number  of  apace  intervals  par  aide. 
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At  Mch  tiM  step  t  and  spae*  point*  ■  on  8^  oq.  (2.7)  ralatM 
p  (t)  to  values  of  p  corraapondlng  to  earlier  tlsei.  Hence  aa  previoitely 
indicated  the  ayatea  of  equation!  obtained  are  aucceaalve  rather  then  aiaul- 
taneous.  Only  at  the  cornera  do  two  adjacent  space  points  have  an  opportuni¬ 
ty  to  Interact  within  a  single  tine  Interval.  Consequently,  at  each  comer,  at 
each  tine  step  there  occur  two  sinultaneous  equations.  The  nunerlcal  results 
obtained  for  the  choice  of  eight  space  steps  on  each  side  of  the  box  and 
sixteen  time  steps  corresponding  to  one  transit  time  are  listed  in  tables  I 
and  II  and  figs.  (3,H  ). 

4.  Discussion  of  Results. 

In  the  case  of  the  box,  the  solution  in  certain  regions  of  space¬ 
time  may  be  identified  with  known  solutions  obtained  by  geometric  acoustics 
[3].  For  example,  for  t  the  solution  along  g  coincides  with  the 

solution  for  an  infinite  wedge  of  vertex  angle  x/2  formed  by  p.  Conse¬ 
quently  the  results  for  S,  »  for  the  first  eight  time  steps  may  be  compared 
with  the  analytic  solution  given  in  (3l>  Since  in  this  case  the  solution 
depends  only  upon  ^  ^  •  comparison  of  values  at  ell  time  steps  may  be 

made  simultaneously  and  the  convergence  of  the  values  for  Increasing  time 
may  be  examined.  The  results  for  the  eighth  tlsw  step  are  plotted  along  with 
the  exact  solution  In  fig.  (  3  )  agreement  Is  excellent j  the 

greatest  error  Is  less  than  3^. 

A  further  estimate  of  the  accuracy  of  the  numerical  solution  may 
be  gained  by  a  comparison  with  the  rigid  box  solution  given  In  ref.  [3]* 

Again  the  agreement  Is  excellent. 

*The  flel^  point  Is  taken  in  the  center  of  the  space  interval  at  /5 
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? 


Bxuilnaeloa  of  fig.  (  H  )  lB<iicot«  that  tha  ataady  atata  aolu* 
tion  of  pal  avary«diara  appaara  to  ba  approaehad  rapidly  avan  within  tha 
intarval  of  ona  transit  tiwa. 
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TABLB  or  VALUES  ON  FRONT  FACI  OP  RIGID  BOX 
UNDER  SYMMETRIC  PULSE:  p(e )  - 

t 

1  1.40 

2  1.36  2.00 


3 

1.34 

1.55 

2.00 

4 

1.32 

1.38 

1.67 

2.00 

5 

1.36 

1.35 

1.50 

1.88 

2.00 

6 

1.34 

1.42 

1.38 

1.64 

2.00 

2.00 

7 

1.33 

1.36 

1.45 

1.44 

1.73 

2.00 

2.00 

8 

1.34 

1.38 

1.39 

1.48 

1.53 

1.83 

2.00 

2.00 

9 

1.34 

1.35 

1.41 

1.44 

1.51 

1.65 

2.00 

1.87 

10 

1.34 

1.36 

1.40 

1.43 

1.47 

1.51 

1.75 

1.62 

11 

1.34 

1.35 

1.38 

1.39 

1.47 

1.52 

1.36 

1.33 

12 

1.34 

1.36 

1.37 

1.40 

1.42 

1.59 

1.24 

1.12 

13 

.  1.34 

1.35 

1.37 

1.37 

1.41 

1.24 

1.15 

1.06 

14 

1.34 

1.35 

1.37 

1.39 

1.25 

1.18 

1.04 

1.09 

15 

1.34 

1.35 

1.36 

1.26 

1.22 

1.05 

1.10 

1.07 

16 

1.34 

1.35 

1.31 

1.21 

1.07 

1.09 

1.05 

1.03 

1 

2 

3 

4 

5 

6 

7 

8 
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TABU  OF  VALUBS  ON  RBAB  SUFACB  OF  KX6XD  BOK 
UNDIR  SMOTRIC  FOUB 

9  0.49 

10  0.93 

11  1.17  0.63 


12 

1.20 

0.89 

0.53 

15 

1.07 

0.94 

0.62 

0.07 

14 

1.04 

1.03 

0 

0.85 

0.46 

15 

1.00 

0.96 

0.92 

0.70 

0.36 

16 

0.99 

0.93 

0.86 

0.79 

0.59 

0.24 

9 

10 

11 

12 

13 

14 
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Appendix  1 


Introduce  the  fundeaentel  eolation  to  the  three  dlaenelonel  wove  equetion 


idlieh  eetlefles 


where  Ra(r— (*>y**)  *  field  point, 

a  aouree  point  and  ^  la  the  Dirac  delta  function  In  one,  three  or  four  dlaen- 
alona  as  required. 

The  secondary  "disturbance”  potential  ^  -  ^^alao  satisfies 


the  wave  equation 


ition.:  f,  =  ^  =  O  «t  t-  O 


with  Initial  conditions; 


Consequently  the  following  equation  la  true 


4Vj>i 


idiere  V. 


Is  a  four  dlwenslonal  volusM  In  space  tiae  and  the  right  hand 


side  Is  4«  O  if  the  point  f>r,tat  Is  Included  In  and  0 

If  It  is  not. 

Two  different  regions  of  Integration  are  eaployed.  They  are  both  bounded 
by  the  initial  time  plane  t^  ■  0,  the  Incident  wave  plane 

obstacle  eurface  extended  Into  space  tlsM  and  a  surface  which  will  Include  tiie 
characteristic  cone  (t-t^)*R  ■  0  In  the  voluoe  of  Integration.  Due  to  the 
singular  bdiavlor  of  the  Integrand,  all  contributions  com  frosi  the  surface 
(t«t^}-R  ■  0  within  the  region  of  Integration.  One  region  of  Intonation  Is 
taken  exterior  to  and  the  other  Interior  to  the  obstacle  surface. 
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Ceasidtr  first,  for  oithlMr  rotion  ,  tho  tora 


«  o 


line*  g  and  Vk.  vanish  for  ovary  valua  of  t  ^  and  and 


arabothOat  t^  «  0. 
o 


The  raaalndar  of  tha  Intagral  la 

'  ,C-**  9^.**  - 

-  j*-:r  su^.  L  s  -  -fv 

by  Graen'a  Thaorem  where  S  la  the  surface  bounding  the  physical  apace 

o 


volume  of  integration  V; 


and  n  is  the  outward  nonaal  to  S  . 
o  0 


Adding  the  two  integrals  corresponding  to  these  two  regions  of  integra¬ 
tion  will  give  an  equation  on  in  terms  of  a  discontinuity  in  <6s 
across  the  obstacle  surface  since  the  normal,  n^,  is  opposite  in  sign  in  the 


two  integrations. 


ft. 

to*Q  S0 


where  (  ]  represents  the  discontinuity  value  and  is  taken  toward  the 

interior  of  the  obstacle  in  physieal  space. 

Since  ^  is  continuous,  the  discontinuity  in  must  be  idmtical 
to  the  discontinuity  in  AS  seen  in  the  text,  the  interior  value  of  ^ 

is  sero  and  the  discontinuity  in  is  thersfors  Just  the  value  immediately 


terior  to  the  obstacle  surface 


\a*k^fe  X\ 

(<P-'-P«)  =  j  Jk.  i  1 -  4^  |m.S 

Vo 
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Beeaui*  pf  tha  bahavlor  of  g,  th«  tlat  intogratlon  can  bo  tokon  thru 
tho  turfoco  InttgrotloB  and  carrlad  out  giving  lq«  (  ). 


Appendix  ll*a 


Tha  InDodlata  neighborhood  of  a  point  on  any  curved  surface  nay  be  eon> 
sidered  plane  aa  an  approximation  which  becomes  exact  as  the  neighborhood 
shrinks  to  the  point.  To  find  the  pressure  field  at  a  point  immediately 
behind  the  wave  front  on  an  obstacle  surface,  it  is  therefore  sufficient  to 
consider  the  interaction  of  the  wave  front  with  an  infinite  plane  obstacle. 
The  result  -  obtained  will  of  course  be  equally  useful  in  the  specific  ex« 
ample  of  the  square  box  for  the  two  planes  which  compose  the  "front"  of  the 


box.  Itie  integral  equation  is  still  valid  but  now  is  applied  only  to 

As  r  approaches  a  point  on  x  approaches  sero.  However, 
mutt  also  approach  xero  except  at  the  indeterminate  point  S  -  0,  l.e.  the 


location  of  the  field  point  on  the  surface.  To  avoid  this  indeteninary 
the  point  R  ■  0  is  removed  from  the  integration  by  distorting  the  surface 


in  the  neighborhood  of  R  a  0  from  a  plane  to  a  semisphere  of  radius  p 
about  R  a  0.  Over  this  surface  *  1.  .  The  integration  is  than  in 
two  parts;  the  semisphere  and  the  remainder  of  the  plane  surface.  The 
letter  integration  is  xero  since  is  xero  and  the  former  integra¬ 


tion  is 
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1!  (  ?  * -f i 

Sine*  In  th«  llait,  th«  t«nu  and  ^  ean-b*  raaoved 

froa  th«  Intagration  and  avaluatad  at  ■  t,  R  ■  0. 

*5{t -X  tcs^-ysmjsj;^.,  ■*•  i 

=  a  C  ■t-yim/s)//. 

f ‘■(r,X)  =  a, 

Appandlx  Il>b 

Bahind  tha  wava  front  on  aach  of  tha  raar  aurfaeaa  of  tha  tquara  box 
axlsta  a  ragion  in  apaca  tine  idilch  la  not  influanead  by  any  of  tha  cornara 
of  tha  box.  Tha  Intagratlona  raqulrad  to  datarailna  tha  praaaura  in  thaaa 
raglona  from  Bq.  (  )  are  than  takan  only  ovar  tha  naighboring  front 

aurfaca  with  valuaa  for  p  on  that  aurfaca  corraaponding  to  tha  infinita 
plana  aolution.  For  axa^^la,  conaidar  a  point  on  aurfaca  i  >  3 

»  a+  a  V  +  ^  ifcJY.  lx.I 


with  Intagratlon  limlta  givan  by 


tih«r«  th«  Halts  of  Intogrstlon  sr«  th«  roots  of 


nils  follows  froa  the  fact  that  ^  ^  aust  be  greater  than  sero  while 

y  Bust  be  less  than  a  in  order  that  the  ends  of  the  front  surface  (l.e. 
'ou  ' 

corners)  will  not  play  any  role  In  the  Integration. 


The  double  Integral  Is  reduced  to  a  single  Integral. 

y#u.  hea  ■ 


o' 

Combining  this  Integral  with  the  term  2L  gives 

^  VL  \  iva  )  _  ? 

^  VsA  L  CXV  Itt'Y#)'’)  i 


•  2  -  2  ■  0 

The  region  In  which  this  solution  will  be  vAld  on  the  rear  surface 
is  t  <t/  /2  from  t  corresponding  to  s-y^  ^  0  and 

X  ^  t/  /2  limiting  points  to  those  behind  the  wave  front. 
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Appendlx  III 


4 

The  expressions  for  Che  obtained  by  Integration  are  given  below 
together  with  the  conditions  under  \^ich  each  expression  la  valid.  These  con* 
dltlons  are  obtained  from  the  limits  In  the  Integral  expressions  for  the  ^ 
and  may  be  Interpreted  In  terms  of  geometric  acoustic  theory  as  representing 
the  diffraction  effects  produced  by  the  corners  of  the  box. 


»  {  »  f  ^If  J. 
rr  J  I  St.  JLlij 


-2lt- 


The  coefflci«nt*  ■“•'  ••ti»fy  th« 

following  symmetry  conditions. 


^  »  u 

*  list  CY-w) 

iWcpm)* 

«« 

iKo)  • 

3  \ 

^  (.y  .)  =  ,<^(f.i-K.>l  (<fY«.J 

y  2ni*  1  ^ 

and  a  similar  set  for  .  }X,  (/w)vhere  X  or  Y  *  — ^»/2  and  p  is 

c  (it'  m  m  c 

the  number  of  steps  in  space  per  side. 

Therefore  only, four  sets  of  coefficients  need  be  computed  for  each 
value  of  m.  llie  integral  expressions  for  the  coefficients  are  simplified 
by  using  a  constant  average  value  in  place  of  over  each  region  of  inte¬ 

gration.  As  Rj  can  never  be  zero  on  any  of  the  surfaces  of  Integration  and 
is  practically  constant  over  all  separate  regions  of  Integration  except  those 
near  the  corners  this  approximation  is  quite  good.  Values  near  the  corners 
must  be  obtained  without  this  approximation. 

,  Ar*m  l-m)  =:■  •  (“i  W  I 


K  ^  li-t 


1 

-  i  Jl'-ll--  lilK-Ut  -v  It. -a^)vw  '  (1=;===^ ) 


^  (<l'-a)C'f  )  •  I  «:>  ,  V*-  u.*-  >C( 

O  ;  .1  Cfc-v)*- 
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Appendlx  IV 


It  is  of  Interest  to  note  that  thvt  pressure  at  the  corner  of 
an  infinite  rigid  wedge  of  arbitrary  angle  o<  may  be  found  exactly  from 
the  integral  equation 

f--  <9^ 


Since 


is  zero  on  both  surfaces  for  a  field  point  at  the  corner 


r  «  a  the  only  contribution  comes  from  the  point  R  b  0  as  in  appendix 

o< 

Ila  although  only  a  portion  ^  of  the  semi-sphere  is  required  to 
exclude  the  corner  point  from  the  surface  of  integration  in  this  case. 


fl-,.,  --  [rr, 


which  is  an  agreement  with  the  known  infinite  rigid  wedge  solution  . 


f] 


r  =o 


4 

3 


o/=  Vz 
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FIG.2  BOX  SHAPED  OBSTACLE 
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Attn:  Laboratory  Division  (1) 

Commanding  Officer 

Frankfotd  Arsenal 

Btidesburg  Station 

Philadelphia  37,  Pennsylvania 

Attn:  Laboratory  Division  (1) 

Office  of  Ordnance  Research 
2127  Myrtle  Drive 
Duke  Station 
Durham,  North  Carolina 
Attn;  Division  of  Engineering 

Sciences  (1) 

NAVY 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn;  Op  37  (1) 

Op  36  (1) 

Op  03EG  (1) 

Commandant,  Marine  Corps 
Headquarters,  U.S.  Marine  Corps 
Washington  25,  D.  C.  (1) 

Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  327  (2) 

Code  348  (1) 

Code  371  (1) 

Code  420  (2) 

Code  423  (2) 

Code  442  (2) 

Code  421  (1) 


Commanding  Officer,  Engineer 
Research  &  Development  Lab. 

Fort  Belvoir,  Virginia 

Attn:  Chief,  Tech.  Intelligence  Br.  (1) 

Director,  Waterways  Experiment  Station 
P.O.  Box  631 
Vicksburg,  Miss. 

Attn:  Library  (1) 

Commanding  Officer 
Engineer  Research  Development 
Laboratory 

Fort  Belvoir,  Virginia  (1) 

Office  of  the  Chief  of  Ordnance 
Department  of  the  Army 
Washington  26,  D.  C. 

Attn:  Research  and  Materials 

Branch  (Ord.  R&D  Div.)  (1) 

ORDTX  -  AR  (1) 


Chief,  Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  AD -2 
AD-22 
RS-7 
TD-42 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washinrton  25,  D.  C. 

Attn:  Ad3 
Re 
Re3 
Re5 
ReN 

Director  of  Naval  Intelligence 
Navy  Department 
Washington  25,  D.  C. 

Attn:  Op-922V 


(1) 

(1) 

(1) 

(1) 


(1) 

(1) 

(1) 

(1) 

(1) 


(1) 


Office  of  the  Chief  Signal  Officer 
Department  of  the  Army 
Washington  25,  D,  C. 

Attn:  Engineering&  Technical  Div.  (1) 


Chief,  Bureau  of  Yards  and  Docks 


Department  of  the  Navy 

Washington  25,  D.  C, 

Attn;  Code  D*213 

(1) 

Code  D-222 

(1) 

Code  D-410C 

(1) 

Code  D-440 

(1) 

Code  D-500 

(1) 
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Commanding  Officer  and  Director 
David  Taylor  Model  Basin 
Washington  7,  D.  C. 

Attn:  Code  140  (1) 

Code  600  (1) 

Code  700  (1) 

Code  720  (1) 

Code  725  (1) 

Code  731  (1) 

Code  740  (2) 

Commander 

U.S.  Naval  Ordnance  Laboratory 
White  Oak,  Maryland 
Atm:  Technical  Library  (2) 

Technical  Evaluation 
Department  (1) 

EE  (1) 

EH  (1) 

R  (1) 

Director 

Material  Laboratory 

New  York  Naval  Shipyard 

Brooklyn  1,  New  York  (1) 

Commanding  Officer  and  Director 
U.S.  Naval  Electronics  Laboratory 
San  Diego  52,  California  (1) 

Attn:  Code  4223  (1) 

Officer-in-Charge 
Naval  Civil  Engineering  Research 
Naval  Civil  Engineering  Research 
and  Evaluation  Laboratory 
U.S.  Naval  Construction  Battalion 


Commanding  Officer  and  Director 
Naval  Engineering  Experiment  Station 
Annapolis,  Maryland  (1) 

Superintendent 

Naval  Post  Graduate  School 

Monterey,  California  (1) 

Commandant 
Marine  Corps  Schools 
Quantico,  Virginia 
Attn:  Director,  Marine  Corps 

Development  Center  (1) 

AIR  FORCE 

Commanding  General 
U.  S.  Ait  Force 
Washington  25,  D.  C. 

Attn:  Research  and  Development 

Division  (1) 

Commander 
Ait  Material  Command 
Wright -Patterson  Air  Force  Base 
Dayton,  Ohio 

Attn:  MCREX-B  (1) 

Structures  Division  (1) 

WCOSI  (1) 

Commander 

U.  S.  Ait  Forece  Institute  of  Technology 
Wti^t -Patterson  Ait  Force  Base 
Dayton,  Ohio 

Attn;  Chief,  Applied  Mechanics 

Group  (1) 


Commandant 
U.S.  Cost  Guard 
1300  E  Street,  N.W. 
Washington  25,  D.  C. 

Attn:  Chief,  Testing  and 

Development  Division 


(1) 


U.S.  Maritime  Administration 
General  Administration  Office  Bldg. 
Washington  25,  D.  C. 

Attn:  Chief,  Division  of  Preliminary 
Design  (1) 

National  Advisory  Committee  for 

Aeronautics 

1512  H  Street,  N.  W. 

Washington  25,  D.  C. 

Attn: 

Director 

Langley  Aeronautical  Laboratory 

Langley  Field,  Virginia 

Attn:  Structures  Division  (2) 

Director 

Forest  Products  Laboratory 
Madison,  Wisconsin  (1) 

Civil  Aeronautics  Administration 
Department  of  Commerce 
Washington  25,  D.  C. 

Attn:  Chief,  Airframe  and  Equipment 
Branch  (1) 

National  Sciences  Foundation 
1520  H  Street,  N.W, 


Port  Hueneme,  California 

(2) 

Director  of  Intelligence 

Attn:  Engineering  Sciences  Div. 

(1) 

Attn;  Code  753 

(1) 

Headquarters,  U.S.  Air  Force 
Washington  25,  D.  C 

National  Academy  of  Science 

Director 

Attn:  P.V.  Branch(Air  TargetsDiv.) 

(1) 

2101  Constitution  Avenue 

Naval  Air  Experimental  Station 

AFOIN-1B2 

(2) 

Washington  25,  D.  C. 

Naval  Ait  Material  Center 

Attn;  Technical  Director,  Committee 

Naval  Base 

Commander 

on  Ships'  Structural  Design 

(1) 

Philadelphia  12,  Pennsylvania 

Air  Research  and  Development  Command 

Executive  Secretary,  Committee 

Attn;  Materials  Laboratory 

(1) 

P.  O.  Box  1395 

on  Undersea  Warfare 

(1) 

Structures  Laboratory 

(1) 

Baltimore  3,  Maryland 

Attn:  RDMPE 

(1) 

Director,  Operations  Research  Office 

Officer-in-Charge 

Johns  Hopkins  University 

Underwater  Explosion  Research 

Commander 

7100  Connecticut  Avenue 

Division 

Wright  Air  Development  Center 

Chevy  Chase,  Maryland 

Norfolk  Naval  Shipyard 

Wright-Patterson  Air  Force  Base 

Washington  15,  D.  C. 

(1) 

Portsmouth,  Virginia 

Dayton,  Ohio 

Attn:  Dr.  A.  H.  Keil 

(2) 

Attn:  Dynamics  Branch 

(1) 

Dr.  Alvin  C.  Graves,  Director 

Aircraft  Laboratory 

(1) 

J-Division,  Los  Alamos  Scientific 

Commander 

WCLSY 

(1) 

Laboratory 

U.S.  Naval  Providing  Grounds 

P.O.  Box  1663 

Dahlgten,  Virginia 

(1) 

Commanding  Officer 

USNNOEU 

Los  Alamos,  New  Mexico 

(1) 

Superintendent 

Kittland  Air  Force  Base 

U.S.  Atomic  Energy  Commission 

Naval  Gun  Factory 

Albuquerque,  New  Mexico 

Classified  Technical  Library 

Washington  25,  D.  C. 

(1) 

Attn:  Code  20 

Technical  Information  Service 

(Dr.  J.N.  Brennan) 

(1) 

1901  Constitution  Avenue,  N.W. 

Commander 

Washington,  D.  C. 

Naval  Ordnance  Test  Station 

PART  C:  OTHER  GOVERNMENT 

Attn:  Mrs.  Jean  M.  O'Leary 

Inyokern,  China  Lake,  California 
Attn:  Physics  Division 

(1) 

ACTIVITIES 

For:  Dr.  Paul  C.  Fine 

(1) 

Mechanics  Branch 

(1) 

U.S.  Atomic  Energy  Commission 

U.S.  Atomic  Energy  Commission 

Washington  25,  D,  C. 

Classified  Technical  Library 

Commander 

Attn;  Director  of  Research 

(2) 

Technical  Information  Service 

Naval  Ordnance  Test  Station 
Underwater  Ordnance  Division 
3202  E.  Foothill  Boulevard 
Pasadena  8,  California 
Attn:  Structures  Division 


(1) 


Director 

National  Bureau  of  Standards 
Washington  25,  D.  C. 

Attn:  Division  of  Mechanics  (1) 

Engineering  Mechanics  Sec.  (1) 
Aircraft  Structures  (1) 


1901  Constitution  Avenue,  N.W. 
Washington,  D.  C, 

Attn:  Mrs.  Jean  M.  O’Leary  (1) 

Sandia  Corporation,  Sandia  Base 

Albuquerque,  New  Mexico 

Attn:  Dr.  Walter  A.  MacNair  (1) 
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Legislative  Reference  Service 
Library  of  Congress 
Washington  25,  D.  C. 

Attn:  Dr.  E.  Wenk  (1) 

U.S.  Atomic  Energy  Commission 
1901  Constitution  Avenue,  N,  W. 
Washington,  D.C.  (1) 

Massachusetts  Institute  of  Technology 

Cambridge  39,  Massachusetts 

Attn:  Dr.  Charles  H,  Norris  (1) 

PART  D:  INVESTIGATORS  ACTIVELY 
ENGAGED  IN  RELATED  RESEARCH 

Professor  Lynn  S.  Beedle 
Fritz  Engineering  Laboratory 
Lehigh  University 

Bethlehem,  Pennsylvania  (1) 

Professor  R.  L.  Bisplinghoff 
Dept,  of  Aeronautical  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts  (1) 

Professor  H.  H.  Bleich 

Department  of  Civil  Engineering 

Columbia  University 

New  York  27.  N.  Y.  (1) 

Professor  B.A.  Boley 

Department  of  Civil  Engineering 

Columbia  University 

New  York,  N.Y.  (1) 

Professor  G.F,  Cartier 
Pierce  Hall 
Harvard  University 

Cambridge  38,  Massachusetts  (1) 

Professor  J.E.  Cermak 
Dept,  of  Civil  Engineering 
Colorado  State  Univ. 

Fort  Collins,  Colorado 

Professor  Herbert  Detesiewicz 
Department  of  Civil  Engineering 
Columbia  University 
632  W.  125th  Street 

New  York  27,  N.  Y.  (1) 

Professor  Lloyd  Donnell 
Department  of  Mechanics 
Illinois  Institute  of  Technology 
Technology  Center 

Chicago  16,  Illinois  (1) 


Professor  D.  C.  Drucker,  Chairman 
Division  of  Engineering 
Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  A.  C.  Eringen 

Department  of  Aeronautical  Engineering 

Purdue  University 

Lafayette,  Indiana  (1) 

Professor  W.  Flugge 

Department  of  Mechanical  Engineering 
Stanford  University  ■ 

Stanford,  California  (1) 


Mr.  Martin  Goland 

Midwest  Research  Institute 

4049  Pennsylvania  Avenue 

Kansas  City  2,  Missouri  (1) 

Professor  J.  N.  Goodier 

Dept,  of  Mechanical  Engineering 

Stanford  University 

Stanford,  California  (1) 

Professor  L.  E.  Goodman 
Engineer  Experiment  Station 
University  of  Minnesota 
Minneapolis,  Minnesota  (1) 

Professor  W.  J.  Hall 
Dept,  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois 

Professor  M.  Hetenyi 

The  Technological  Institute 

Northwestern  University 

Evanston,  Illinois  (1) 

Professor  P.  G.  Hodge 

Department  of  Aeronautical  Engineering 

and  Applied  Mechanics 

Polytechnic  Institute  of  Brooklyn 

99  Livingston  Street 

Brooklyn  2,  N.  Y. 

Professor  N.  J.  Hoff 
Stanford  University 

Stanford,  California  (1) 

Professor  W.  H.  Hoppmann,  II 
Department  of  Mechanical  Engineering 
Johns  Hopkins  University 
Baltimore,  Maryland  (1) 

Professor  Bruce  G.  Johnston 

University  of  Michigan 

Ann  Arbor,  Michigan  (1) 

Professor  J.  Kempner 

Department  of  Aeronautical  Engineering 

and  Applied  Mechanics 

Polytechnic  Institute  of  Brooklyn 

99  Livingston  Street 

Brooklyn  2,  New  York  (1) 

Professor  H.  L.  Langhaar 

Department  of  Theoretical  and 

Applied  Mechanics 

University  of  Illinois 

Urbana,  Illinois  (1) 

Professor  B,J.  Lazan,  Director 
Engineering  Experiment  Station 
University  of  Minnesota 
Minneapolis  14,  Minnesota  (1) 

Professor  E.  H.  Lee 

Division  of  Applied  Mathematics 

Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  George  H.  Lee 
Director  of  Research 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  (1) 


Mr.  M.  M.  Lemcoe 
Southwest  Research  Institute 
8500  Culebra  Road 

San  Antonio  6,  Texas  (1) 

Professor  Paul  Lieber 
Geology  Department 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  (1) 

Professor  Hsu  Lo 
School  of  Engineering 
Purdue  University 

Lafayette,  Indiana  (1) 

Professor  R.  D.  Mindlin 

Department  of  Civil  Engineering 

Columbia  University 

632  W.  T25th  Street 

New  York  27,  N.  Y.  (1) 

Dr.  A.  Nadai 

136  Cherry  Valley  Road 

Pittsburgh  21,  Pennsylvania  (1) 

Professor  Paul  M.  Naghdi 
Department  of  Mechanical 
Engineering 

University  of  California 

Berkeley,  California  (1) 

Professor  William  A.  Nash 
Department  of  Engineering 
Mechanics 
University  of  Florida 
Gainesville,  Florida  (1) 

Professor  N.  M.  Newmark 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  Aris  Phillips 
Department  of  Civil  Engineering 
15  Prospect  Street 
Yale  University 

New  Haven,  Connecticut  (1) 

Professor  W.  Prager,  Chairman 
Physical  Sciences  Council 
Brown  University 

Providence  12,  Rhode  Island  (1) 

Professor  E.  Reissner 
Department  of  Mathematics 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts  (1) 

Professor  M.  A.  Sadowsky 
Department  of  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York 

Dr.  B.  W.  Shaffer 

Dept,  of  Mechanical  Engineering 

New  York  University 

45  Fourth  Avenue 

New  York  53,  New  York  (1) 

Professor  C.  B.  Smith 

College  of  Arts  and  Sciences 

Department  of  Mathematics 

Walker  Hall 

University  of  Florida 

Gainesville,  Florida  (1) 
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Professor  J.  Stallmeyer 
Department  of  Civil  Engineering 
University  of  Illinois 
Urban,  Illinois 

Professor  Eli  Sternberg 
Brown  University 

Providence,  Rhode  Island  (1) 

Professor  S.P.  Timoshenko 
School  of  Engineering 
Stanford  University 

Stanford,  California  (1) 

Professor  A. A,  Velestos 
Department  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois  (1) 

Professor  Enrico  Voltetra 
Department  of  Engineering  Mech 
Mechanics 
University  of  Texas 
Austin  T2,  Texas  (1) 

Professor  Dana  Young 
Yale  University 

New  Haven,  Connecticut  (1) 

Project  Staff  (10) 

For  your  future  distribution  (10) 
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